Disentangling the relationships of body mass index and circulating sex hormone concentrations in mammographic density using Mendelian randomization.

PURPOSE: Mammographic density phenotypes, adjusted for age and body mass index (BMI), are strong predictors of breast cancer risk. BMI is associated with mammographic density measures, but the role of circulating sex hormone concentrations is less clear. We investigated the relationship between BMI, circulating sex hormone concentrations, and mammographic density phenotypes using Mendelian randomization (MR). METHODS: We applied two-sample MR approaches to assess the association between genetically predicted circulating concentrations of sex hormones [estradiol, testosterone, sex hormone-binding globulin (SHBG)], BMI, and mammographic density phenotypes (dense and non-dense area). We created instrumental variables from large European ancestry-based genome-wide association studies and applied estimates to mammographic density phenotypes in up to 14,000 women of European ancestry. We performed analyses overall and by menopausal status. RESULTS: Genetically predicted BMI was positively associated with non-dense area (IVW: ß = 1.79; 95% CI = 1.58, 2.00; p = 9.57 × 10-63) and inversely associated with dense area (IVW: ß = - 0.37; 95% CI = - 0.51,- 0.23; p = 4.7 × 10-7). We observed weak evidence for an association of circulating sex hormone concentrations with mammographic density phenotypes, specifically inverse associations between genetically predicted testosterone concentration and dense area (ß = - 0.22; 95% CI = - 0.38, - 0.053; p = 0.009) and between genetically predicted estradiol concentration and non-dense area (ß = - 3.32; 95% CI = - 5.83, - 0.82; p = 0.009), although results were not consistent across a range of MR approaches. CONCLUSION: Our findings support a positive causal association between BMI and mammographic non-dense area and an inverse association between BMI and dense area. Evidence was weaker and inconsistent for a causal effect of circulating sex hormone concentrations on mammographic density phenotypes. Based on our findings, associations between circulating sex hormone concentrations and mammographic density phenotypes are weak at best.

Temporal changes in mammographic breast density and breast cancer risk among women with benign breast disease.

INTRODUCTION: Benign breast disease (BBD) and high mammographic breast density (MBD) are prevalent and independent risk factors for invasive breast cancer. It has been suggested that temporal changes in MBD may impact future invasive breast cancer risk, but this has not been studied among women with BBD. METHODS: We undertook a nested case-control study within a cohort of 15,395 women with BBD in Kaiser Permanente Northwest (KPNW; 1970-2012, followed through mid-2015). Cases (n = 261) developed invasive breast cancer > 1 year after BBD diagnosis, whereas controls (n = 249) did not have breast cancer by the case diagnosis date. Cases and controls were individually matched on BBD diagnosis age and plan membership duration. Standardized %MBD change (per 2 years), categorized as stable/any increase (≥ 0%), minimal decrease of less than 5% or a decrease greater than or equal to 5%, was determined from baseline and follow-up mammograms. Associations between MBD change and breast cancer risk were examined using adjusted unconditional logistic regression. RESULTS: Overall, 64.5% (n = 329) of BBD patients had non-proliferative and 35.5% (n = 181) had proliferative disease with/without atypia. Women with an MBD decrease (≤ - 5%) were less likely to develop breast cancer (Odds Ratio (OR) 0.64; 95% Confidence Interval (CI) 0.38, 1.07) compared with women with minimal decreases. Associations were stronger among women ≥ 50 years at BBD diagnosis (OR 0.48; 95% CI 0.25, 0.92) and with proliferative BBD (OR 0.32; 95% CI 0.11, 0.99). DISCUSSION: Assessment of temporal MBD changes may inform risk monitoring among women with BBD, and strategies to actively reduce MBD may help decrease future breast cancer risk.

Host, reproductive, and lifestyle factors in relation to quantitative histologic metrics of the normal breast.

BACKGROUND: Emerging data indicate that variations in quantitative epithelial and stromal tissue composition and their relative abundance in benign breast biopsies independently impact risk of future invasive breast cancer. To gain further insights into breast cancer etiopathogenesis, we investigated associations between epidemiological factors and quantitative tissue composition metrics of the normal breast. METHODS: The study participants were 4108 healthy women ages 18-75 years who voluntarily donated breast tissue to the US-based Susan G. Komen Tissue Bank (KTB; 2008-2019). Using high-accuracy machine learning algorithms, we quantified the percentage of epithelial, stromal, adipose, and fibroglandular tissue, as well as the proportion of fibroglandular tissue that is epithelium relative to stroma (i.e., epithelium-to-stroma proportion, ESP) on digitized hematoxylin and eosin (H&E)-stained normal breast biopsy specimens. Data on epidemiological factors were obtained from participants using a detailed questionnaire administered at the time of tissue donation. Associations between epidemiological factors and square root transformed tissue metrics were investigated using multivariable linear regression models. RESULTS: With increasing age, the amount of stromal, epithelial, and fibroglandular tissue declined and adipose tissue increased, while that of ESP demonstrated a bimodal pattern. Several epidemiological factors were associated with individual tissue composition metrics, impacting ESP as a result. Compared with premenopausal women, postmenopausal women had lower ESP [ß (95% Confidence Interval (CI)) = -0.28 (- 0.43, - 0.13); P < 0.001] with ESP peaks at 30-40 years and 60-70 years among pre- and postmenopausal women, respectively. Pregnancy [ß (95%CI) vs nulligravid = 0.19 (0.08, 0.30); P < 0.001] and increasing number of live births (P-trend < 0.001) were positively associated with ESP, while breastfeeding was inversely associated with ESP [ß (95%CI) vs no breastfeeding = -0.15 (- 0.29, - 0.01); P = 0.036]. A positive family history of breast cancer (FHBC) [ß (95%CI) vs no FHBC = 0.14 (0.02-0.26); P = 0.02], being overweight or obese [ß (95%CI) vs normal weight = 0.18 (0.06-0.30); P = 0.004 and 0.32 (0.21-0.44); P < 0.001, respectively], and Black race [ß (95%CI) vs White = 0.12 (- 0.005, 0.25); P = 0.06] were positively associated with ESP. CONCLUSION: Our findings revealed that cumulative exposure to etiological factors over the lifespan impacts normal breast tissue composition metrics, individually or jointly, to alter their dynamic equilibrium, with potential implications for breast cancer susceptibility and tumor etiologic heterogeneity.

Outdoor air pollution and histologic composition of normal breast tissue.

BACKGROUND: Biologic pathways underlying the association between outdoor air pollution and breast cancer risk are poorly understood. Breast tissue composition may reflect cumulative exposure to breast cancer risk factors and has been associated with breast cancer risk among patients with benign breast disease. Herein, we evaluated whether fine particulate matter (PM2.5) was associated with the histologic composition of normal breast tissue. METHODS: Machine-learning algorithms were applied to digitized hematoxylin and eosin-stained biopsies of normal breast tissue to quantify the epithelium, stroma, adipose and total tissue area from 3,977 individuals aged 18-75 years from a primarily Midwestern United States population who donated breast tissue samples to the Susan G. Komen Tissue Bank (2009-2019). Annual levels of PM2.5 were assigned to each woman's residential address based on year of tissue donation. We applied predictive k-means to assign participants to clusters with similar PM2.5 chemical composition and used linear regression to examine the cross-sectional associations between a 5-µg/m3 increase in PM2.5 and square root-transformed proportions of epithelium, stroma, adipose, and epithelium-to-stroma proportion [ESP], overall and by PM2.5 cluster. RESULTS: Higher residential PM2.5 was associated with lower proportion of breast stromal tissue [ß = -0.93, 95% confidence interval: (-1.52, -0.33)], but was not related to the proportion of epithelium [ß = -0.11 (-0.34, 0.11)]. Although PM2.5 was not associated with ESP overall [ß = 0.24 (-0.16, 0.64)], the association significantly differed by PM2.5 chemical composition (p-interaction = 0.04), with a positive association evident only among an urban, Midwestern cluster with higher concentrations of nitrate (NO3-) and ammonium (NH4+) [ß = 0.49 (0.03, 0.95)]. CONCLUSIONS: Our findings are consistent with a possible role of PM2.5 in breast cancer etiology and suggest that changes in breast tissue composition may be a potential pathway by which outdoor air pollution impacts breast cancer risk. This study further underscores the importance of considering heterogeneity in PM2.5 composition and its impact on breast carcinogenesis.

Genome-wide and transcriptome-wide association studies of mammographic density phenotypes reveal novel loci.

BACKGROUND: Mammographic density (MD) phenotypes, including percent density (PMD), area of dense tissue (DA), and area of non-dense tissue (NDA), are associated with breast cancer risk. Twin studies suggest that MD phenotypes are highly heritable. However, only a small proportion of their variance is explained by identified genetic variants. METHODS: We conducted a genome-wide association study, as well as a transcriptome-wide association study (TWAS), of age- and BMI-adjusted DA, NDA, and PMD in up to 27,900 European-ancestry women from the MODE/BCAC consortia. RESULTS: We identified 28 genome-wide significant loci for MD phenotypes, including nine novel signals (5q11.2, 5q14.1, 5q31.1, 5q33.3, 5q35.1, 7p11.2, 8q24.13, 12p11.2, 16q12.2). Further, 45% of all known breast cancer SNPs were associated with at least one MD phenotype at p < 0.05. TWAS further identified two novel genes (SHOX2 and CRISPLD2) whose genetically predicted expression was significantly associated with MD phenotypes. CONCLUSIONS: Our findings provided novel insight into the genetic background of MD phenotypes, and further demonstrated their shared genetic basis with breast cancer.

Rapid Reductions in Breast Density following Tamoxifen Therapy as Evaluated by Whole-Breast Ultrasound Tomography.

PURPOSE: Women whose mammographic breast density declines within 12-18 months of initiating tamoxifen for chemoprevention or adjuvant treatment show improved therapeutic responses compared with those whose density is unchanged. We tested whether measuring changes in sound speed (a surrogate of breast density) using ultrasound tomography (UST) could enable rapid identification of favorable responses to tamoxifen. METHODS: We evaluated serial density measures at baseline and at 1 to 3, 4 to 6, and 12+ months among 74 women (aged 30-70 years) following initiation of tamoxifen for clinical indications, including an elevated risk of breast cancer (20%) and diagnoses of in situ (39%) or invasive (40%) breast carcinoma, enrolled at Karmanos Cancer Institute and Henry Ford Health System (Detroit, MI, USA). For comparison, we evaluated an untreated group with screen negative mammography and frequency-matched on age, race, and menopausal status (n = 150), at baseline and 12 months. Paired t-tests were used to assess differences in UST sound speed over time and between tamoxifen-treated and untreated patients. RESULTS: Sound speed declined steadily over the 12 month period among patients receiving tamoxifen (mean (SD): -3.0 (8.2) m/s; p = 0.001), whereas density remained unchanged in the untreated group (mean (SD): 0.4 (7.1) m/s; p = 0.75 (relative change between groups: p = 0.0009)). In the tamoxifen group, we observed significant sound speed reductions as early as 4-6 months after tamoxifen initiation (mean (SD): -2.1 (6.8) m/s; p = 0.008). Sound speed reductions were greatest among premenopausal patients (P-interaction = 0.0002) and those in the middle and upper tertiles of baseline sound speed (P-interaction = 0.002). CONCLUSIONS: UST can image rapid declines in sound speed following initiation of tamoxifen. Given that sound speed and mammographic density are correlated, we propose that UST breast imaging may capture early responses to tamoxifen, which in turn may have utility in predicting therapeutic efficacy.

Mammary collagen architecture and its association with mammographic density and lesion severity among women undergoing image-guided breast biopsy.

BACKGROUND: Elevated mammographic breast density is a strong breast cancer risk factor with poorly understood etiology. Increased deposition of collagen, one of the main fibrous proteins present in breast stroma, has been associated with increased mammographic density. Collagen fiber architecture has been linked to poor outcomes in breast cancer. However, relationships of quantitative collagen fiber features assessed in diagnostic biopsies with mammographic density and lesion severity are not well-established. METHODS: Clinically indicated breast biopsies from 65 in situ or invasive breast cancer cases and 73 frequency matched-controls with a benign biopsy result were used to measure collagen fiber features (length, straightness, width, alignment, orientation and density (fibers/µm2)) using second harmonic generation microscopy in up to three regions of interest (ROIs) per biopsy: normal, benign breast disease, and cancer. Local and global mammographic density volumes were quantified in the ipsilateral breast in pre-biopsy full-field digital mammograms. Associations of fibrillar collagen features with mammographic density and severity of biopsy diagnosis were evaluated using generalized estimating equation models with an independent correlation structure to account for multiple ROIs within each biopsy section. RESULTS: Collagen fiber density was positively associated with the proportion of stroma on the biopsy slide (p < 0.001) and with local percent mammographic density volume at both the biopsy target (p = 0.035) and within a 2 mm perilesional ring (p = 0.02), but not with global mammographic density measures. As severity of the breast biopsy diagnosis increased at the ROI level, collagen fibers tended to be less dense, shorter, straighter, thinner, and more aligned with one another (p < 0.05). CONCLUSIONS: Collagen fiber density was positively associated with local, but not global, mammographic density, suggesting that collagen microarchitecture may not translate into macroscopic mammographic features. However, collagen fiber features may be markers of cancer risk and/or progression among women referred for biopsy based on abnormal breast imaging.

Breast cancer risk factors in relation to molecular subtypes in breast cancer patients from Kenya.

BACKGROUND: Few studies have investigated risk factor heterogeneity by molecular subtypes in indigenous African populations where prevalence of traditional breast cancer (BC) risk factors, genetic background, and environmental exposures show marked differences compared to European ancestry populations. METHODS: We conducted a case-only analysis of 838 pathologically confirmed BC cases recruited from 5 groups of public, faith-based, and private institutions across Kenya between March 2012 to May 2015. Centralized pathology review and immunohistochemistry (IHC) for key markers (ER, PR, HER2, EGFR, CK5-6, and Ki67) was performed to define subtypes. Risk factor data was collected at time of diagnosis through a questionnaire. Multivariable polytomous logistic regression models were used to determine associations between BC risk factors and tumor molecular subtypes, adjusted for clinical characteristics and risk factors. RESULTS: The median age at menarche and first pregnancy were 14 and 21 years, median number of children was 3, and breastfeeding duration was 62 months per child. Distribution of molecular subtypes for luminal A, luminal B, HER2-enriched, and triple negative (TN) breast cancers was 34.8%, 35.8%, 10.7%, and 18.6%, respectively. After adjusting for covariates, compared to patients with ER-positive tumors, ER-negative patients were more likely to have higher parity (OR = 2.03, 95% CI = (1.11, 3.72), p = 0.021, comparing ≥ 5 to ≤ 2 children). Compared to patients with luminal A tumors, luminal B patients were more likely to have lower parity (OR = 0.45, 95% CI = 0.23, 0.87, p = 0.018, comparing ≥ 5 to ≤ 2 children); HER2-enriched patients were less likely to be obese (OR = 0.36, 95% CI = 0.16, 0.81, p = 0.013) or older age at menopause (OR = 0.38, 95% CI = 0.15, 0.997, p = 0.049). Body mass index (BMI), either overall or by menopausal status, did not vary significantly by ER status. Overall, cumulative or average breastfeeding duration did not vary significantly across subtypes. CONCLUSIONS: In Kenya, we found associations between parity-related risk factors and ER status consistent with observations in European ancestry populations, but differing associations with BMI and breastfeeding. Inclusion of diverse populations in cancer etiology studies is needed to develop population and subtype-specific risk prediction/prevention strategies.

Relation of Quantitative Histologic and Radiologic Breast Tissue Composition Metrics With Invasive Breast Cancer Risk.

Background: Benign breast disease (BBD) is a strong breast cancer risk factor, but identifying patients that might develop invasive breast cancer remains a challenge. Methods: By applying machine-learning to digitized hematoxylin and eosin-stained biopsies and computer-assisted thresholding to mammograms obtained circa BBD diagnosis, we generated quantitative tissue composition metrics and determined their association with future invasive breast cancer diagnosis. Archival breast biopsies and mammograms were obtained for women (18-86 years of age) in a case-control study, nested within a cohort of 15 395 BBD patients from Kaiser Permanente Northwest (1970-2012), followed through mid-2015. Patients who developed incident invasive breast cancer (ie, cases; n = 514) and those who did not (ie, controls; n = 514) were matched on BBD diagnosis age and plan membership duration. All statistical tests were 2-sided. Results: Increasing epithelial area on the BBD biopsy was associated with increasing breast cancer risk (odds ratio [OR]Q4 vs Q1 = 1.85, 95% confidence interval [CI] = 1.13 to 3.04; P trend = .02). Conversely, increasing stroma was associated with decreased risk in nonproliferative, but not proliferative, BBD (P heterogeneity = .002). Increasing epithelium-to-stroma proportion (ORQ4 vs Q1 = 2.06, 95% CI =1.28 to 3.33; P trend = .002) and percent mammographic density (MBD) (ORQ4 vs Q1 = 2.20, 95% CI = 1.20 to 4.03; P trend = .01) were independently and strongly predictive of increased breast cancer risk. In combination, women with high epithelium-to-stroma proportion and high MBD had substantially higher risk than those with low epithelium-to-stroma proportion and low MBD (OR = 2.27, 95% CI = 1.27 to 4.06; P trend = .005), particularly among women with nonproliferative (P trend = .01) vs proliferative (P trend = .33) BBD. Conclusion: Among BBD patients, increasing epithelium-to-stroma proportion on BBD biopsies and percent MBD at BBD diagnosis were independently and jointly associated with increasing breast cancer risk. These findings were particularly striking for women with nonproliferative disease (comprising approximately 70% of all BBD patients), for whom relevant predictive biomarkers are lacking.

Risk factors for breast cancer development by tumor characteristics among women with benign breast disease.

BACKGROUND: Among women diagnosed with invasive breast cancer, 30% have a prior diagnosis of benign breast disease (BBD). Thus, it is important to identify factors among BBD patients that elevate invasive cancer risk. In the general population, risk factors differ in their associations by clinical pathologic features; however, whether women with BBD show etiologic heterogeneity in the types of breast cancers they develop remains unknown. METHODS: Using a nested case-control study of BBD and breast cancer risk conducted in a community healthcare plan (Kaiser Permanente Northwest), we assessed relationships of histologic features in BBD biopsies and patient characteristics with subsequent breast cancer risk and tested for heterogeneity of associations by estrogen receptor (ER) status, tumor grade, and size. The study included 514 invasive breast cancer cases (median follow-up of 9 years post-BBD diagnosis) and 514 matched controls, diagnosed with proliferative or non-proliferative BBD between 1971 and 2006, with follow-up through mid-2015. Odds ratios (ORs) and 95% confidence intervals (CIs) were obtained using multivariable polytomous logistic regression models. RESULTS: Breast cancers were predominantly ER-positive (86%), well or moderately differentiated (73%), small (74% < 20 mm), and stage I/II (91%). Compared to patients with non-proliferative BBD, proliferative BBD with atypia conferred increased risk for ER-positive cancer (OR = 5.48, 95% CI = 2.14-14.01) with only one ER-negative case, P-heterogeneity = 0.45. The presence of columnar cell lesions (CCLs) at BBD diagnosis was associated with a 1.5-fold increase in the risk of both ER-positive and ER-negative tumors, with a 2-fold increase (95% CI = 1.21-3.58) observed among postmenopausal women (56%), independent of proliferative BBD status with and without atypia. We did not identify statistically significant differences in risk factor associations by tumor grade or size. CONCLUSION: Most tumors that developed after a BBD diagnosis in this cohort were highly treatable low-stage ER-positive tumors. CCL in BBD biopsies may be associated with moderately increased risk, independent of BBD histology, and irrespective of ER status.

Outdoor air pollution and terminal duct lobular involution of the normal breast.

BACKGROUND: Exposure to certain outdoor air pollutants may be associated with a higher risk of breast cancer, though potential underlying mechanisms are poorly understood. We examined whether outdoor air pollution was associated with involution of terminal duct lobular units (TDLUs), the histologic site where most cancers arise and an intermediate marker of breast cancer risk. METHODS: Pathologist-enumerated TDLUs were assessed in H&E (hematoxylin and eosin)-stained breast tissue sections from 1904 US women ages 18-75 who donated to the Susan G. Komen Tissue Bank (2009-2012). The 2009 annual fine particulate matter < 2.5 µm in diameter (PM2.5) total mass (µg/m3) at each woman's residential address was estimated from the Environmental Protection Agency's Downscaler Model combining Community Multiscale Air Quality (CMAQ) System modeling with air quality monitoring data. We secondarily considered CMAQ-modeled components of PM2.5 and gaseous pollutants. We used K-means clustering to identify groups of individuals with similar levels of PM2.5 components, selecting groups via cluster stability analysis. Relative rates (RRs) and 95% confidence intervals (95% CIs) for the association between air pollutants and TDLU counts were estimated from a zero-inflated negative binomial regression model adjusted for potential confounders. RESULTS: PM2.5 total mass was associated with higher TDLU counts among all women (interquartile range (IQR) increase, RR = 1.06; 95% CI: 1.01-1.11). This association was evident among both premenopausal and postmenopausal women (premenopausal RR = 1.05, 95% CI: 1.00-1.11; postmenopausal RR = 1.11, 95% CI: 1.00-1.23). We identified 3 groups corresponding to clusters that varied geographically and roughly represented high, medium, and low levels of PM2.5 components relative to population mean levels. Compared to the cluster with low levels, the clusters with both high (RR = 1.74; 95% CI: 1.08-2.80) and medium (RR = 1.82; 95% CI: 1.13-2.93) levels were associated with higher TDLU counts; although not significantly different, the magnitude of the associations was stronger among postmenopausal women. CONCLUSIONS: Higher PM2.5 levels were associated with reduced TDLU involution as measured by TDLU counts. Air pollution exposure may influence the histologic characteristics of normal tissue which could in turn affect breast cancer risk.

Polygenic risk score for the prediction of breast cancer is related to lesser terminal duct lobular unit involution of the breast.

Terminal duct lobular units (TDLUs) are the predominant anatomical structures where breast cancers originate. Having lesser degrees of age-related TDLU involution, measured as higher TDLUs counts or more epithelial TDLU substructures (acini), is related to increased breast cancer risk among women with benign breast disease (BBD). We evaluated whether a recently developed polygenic risk score (PRS) based on 313-common variants for breast cancer prediction is related to TDLU involution in the background, normal breast tissue, as this could provide mechanistic clues on the genetic predisposition to breast cancer. Among 1398 women without breast cancer, higher values of the PRS were significantly associated with higher TDLU counts (P = 0.004), but not with acini counts (P = 0.808), in histologically normal tissue samples from donors and diagnostic BBD biopsies. Mediation analysis indicated that TDLU counts may explain a modest proportion (≤10%) of the association of the 313-variant PRS with breast cancer risk. These findings suggest that TDLU involution might be an intermediate step in the association between common genetic variation and breast cancer risk.

Using Whole Breast Ultrasound Tomography to Improve Breast Cancer Risk Assessment: A Novel Risk Factor Based on the Quantitative Tissue Property of Sound Speed.

Mammographic percent density (MPD) is an independent risk factor for developing breast cancer, but its inclusion in clinical risk models provides only modest improvements in individualized risk prediction, and MPD is not typically assessed in younger women because of ionizing radiation concerns. Previous studies have shown that tissue sound speed, derived from whole breast ultrasound tomography (UST), a non-ionizing modality, is a potential surrogate marker of breast density, but prior to this study, sound speed has not been directly linked to breast cancer risk. To that end, we explored the relation of sound speed and MPD with breast cancer risk in a case-control study, including 61 cases with recent breast cancer diagnoses and a comparison group of 165 women, frequency matched to cases on age, race, and menopausal status, and with a recent negative mammogram and no personal history of breast cancer. Multivariable odds ratios (ORs) and 95% confidence intervals (CIs) were estimated for the relation of quartiles of MPD and sound speed with breast cancer risk adjusted for matching factors. Elevated MPD was associated with increased breast cancer risk, although the trend did not reach statistical significance (OR per quartile = 1.27, 95% CI: 0.95, 1.70; ptrend = 0.10). In contrast, elevated sound speed was significantly associated with breast cancer risk in a dose-response fashion (OR per quartile = 1.83, 95% CI: 1.32, 2.54; ptrend = 0.0003). The OR trend for sound speed was statistically significantly different from that observed for MPD (p = 0.005). These findings suggest that whole breast sound speed may be more strongly associated with breast cancer risk than MPD and offer future opportunities for refining the magnitude and precision of risk associations in larger, population-based studies, including women younger than usual screening ages.

Relationship of Serum Progesterone and Progesterone Metabolites with Mammographic Breast Density and Terminal Ductal Lobular Unit Involution among Women Undergoing Diagnostic Breast Biopsy.

The association of progesterone/progesterone metabolites with elevated mammographic breast density (MBD) and delayed age-related terminal duct lobular unit (TDLU) involution, strong breast cancer risk factors, has received limited attention. Using a reliable liquid chromatography-tandem mass-spectrometry assay, we quantified serum progesterone/progesterone metabolites and explored cross-sectional relationships with MBD and TDLU involution among women, ages 40-65, undergoing diagnostic breast biopsy. Quantitative MBD measures were estimated in pre-biopsy digital mammograms. TDLU involution was quantified in diagnostic biopsies. Adjusted partial correlations and trends across MBD/TDLU categories were calculated. Pregnenolone was positively associated with percent MBD-area (MBD-A, rho: 0.30; p-trend = 0.01) among premenopausal luteal phase women. Progesterone tended to be positively associated with percent MBD-A among luteal phase (rho: 0.26; p-trend = 0.07) and postmenopausal (rho: 0.17; p-trend = 0.04) women. Consistent with experimental data, implicating an elevated 5α-pregnanes/3α-dihydroprogesterone (5αP/3αHP) metabolite ratio in breast cancer, higher 5αP/3αHP was associated with elevated percent MBD-A among luteal phase (rho: 0.29; p-trend = 0.08), but not postmenopausal women. This exploratory analysis provided some evidence that endogenous progesterone and progesterone metabolites might be correlated with MBD, a strong breast cancer risk factor, in both pre- and postmenopausal women undergoing breast biopsy. Additional studies are needed to understand the role of progesterone/progesterone metabolites in breast tissue composition and breast cancer risk.

Association of lifestyle and clinical characteristics with receipt of radiotherapy treatment among women diagnosed with DCIS in the NIH-AARP Diet and Health Study.

PURPOSE: The long-term risks and benefits of radiotherapy for ductal carcinoma in situ (DCIS) remain unclear. Recent data from the Surveillance, Epidemiology and End Results (SEER) registries showed that DCIS-associated radiotherapy treatment significantly increased risk of second non-breast cancers including lung cancer. To help understand those observations and whether breast cancer risk factors are related to radiotherapy treatment decision-making, we examined associations between lifestyle and clinical factors with DCIS radiotherapy receipt. METHODS: Among 1628 participants from the NIH-AARP Diet and Health Study, diagnosed with incident DCIS (1995-2011), we examined associations between lifestyle and clinical factors with radiotherapy receipt. Radiotherapy and clinical information were ascertained from state cancer registries. Odds ratios (ORs) and 95% confidence intervals (CIs) for radiotherapy receipt (yes/no) were estimated from multivariable logistic regression. RESULTS: Overall, 45% (n = 730) received radiotherapy. No relationships were observed for most lifestyle factors and radiotherapy receipt, including current smoking (OR 0.97, 95%CI 0.70, 1.34). However positive associations were observed for moderate alcohol consumption and infrequent physical activity. The strongest associations were observed for radiotherapy receipt and more recent diagnoses (2005-2011 vs. 1995-1999; OR 1.60, 95%CI 1.14, 2.25), poorly versus well-differentiated tumors (OR 1.69, 95%CI 1.16, 2.46) and endocrine therapy (OR 3.37, 95%CI 2.56, 4.44). CONCLUSIONS: Clinical characteristics were the strongest determinants of DCIS radiotherapy. Receipt was largely unrelated to lifestyle factors suggesting that the previously observed associations in SEER were likely not confounded by these lifestyle factors. Further studies are needed to understand mechanisms driving radiotherapy-associated second malignancies following DCIS, to identify prevention opportunities for this growing population.

Involution of Breast Lobules, Mammographic Breast Density and Prognosis Among Tamoxifen-Treated Estrogen Receptor-Positive Breast Cancer Patients.

Mammographic breast density (MD) reflects breast fibroglandular content. Its decline following adjuvant tamoxifen treated, estrogen receptor (ER)-positive breast cancer has been associated with improved outcomes. Breast cancers arise from structures termed lobules, and lower MD is associated with increased age-related lobule involution. We assessed whether pre-treatment involution influenced associations between MD decline and risk of breast cancer-specific death. ER-positive tamoxifen treated patients diagnosed at Kaiser Permanente Northwest (1990-2008) were defined as cases who died of breast cancer (n = 54) and matched controls (remained alive over similar follow-up; n = 180). Lobule involution was assessed by examining terminal duct lobular units (TDLUs) in benign tissues surrounding cancers as TDLU count/mm2, median span and acini count/TDLU. MD (%) was measured in the unaffected breast at baseline (median 6-months before) and follow-up (median 12-months after tamoxifen initiation). TDLU measures and baseline MD were positively associated among controls (p < 0.05). In multivariable regression models, MD decline (≥10%) was associated with reduced risk of breast cancer-specific death before (odds ratio (OR): 0.41, 95% CI: 0.18-0.92) and after (OR: 0.41, 95% CI: 0.18-0.94) adjustment for TDLU count/mm2, TDLU span (OR: 0.34, 95% CI: 0.14-0.84), and acini count/TDLU (OR: 0.33, 95% CI: 0.13-0.81). MD decline following adjuvant tamoxifen is associated with reduced risk of breast cancer-specific death, irrespective of pre-treatment lobule involution.

Application of convolutional neural networks to breast biopsies to delineate tissue correlates of mammographic breast density.

Breast density, a breast cancer risk factor, is a radiologic feature that reflects fibroglandular tissue content relative to breast area or volume. Its histology is incompletely characterized. Here we use deep learning approaches to identify histologic correlates in radiologically-guided biopsies that may underlie breast density and distinguish cancer among women with elevated and low density. We evaluated hematoxylin and eosin (H&E)-stained digitized images from image-guided breast biopsies (n = 852 patients). Breast density was assessed as global and localized fibroglandular volume (%). A convolutional neural network characterized H&E composition. In total 37 features were extracted from the network output, describing tissue quantities and morphological structure. A random forest regression model was trained to identify correlates most predictive of fibroglandular volume (n = 588). Correlations between predicted and radiologically quantified fibroglandular volume were assessed in 264 independent patients. A second random forest classifier was trained to predict diagnosis (invasive vs. benign); performance was assessed using area under receiver-operating characteristics curves (AUC). Using extracted features, regression models predicted global (r = 0.94) and localized (r = 0.93) fibroglandular volume, with fat and non-fatty stromal content representing the strongest correlates, followed by epithelial organization rather than quantity. For predicting cancer among high and low fibroglandular volume, the classifier achieved AUCs of 0.92 and 0.84, respectively, with epithelial organizational features ranking most important. These results suggest non-fatty stroma, fat tissue quantities and epithelial region organization predict fibroglandular volume. The model holds promise for identifying histological correlates of cancer risk in patients with high and low density and warrants further evaluation.

Using Digital Pathology to Understand Epithelial Characteristics of Benign Breast Disease among Women Undergoing Diagnostic Image-Guided Breast Biopsy.

Delayed terminal duct lobular unit (TDLU) involution is associated with elevated mammographic breast density (MD). Both are independent breast cancer risk factors among women with benign breast disease (BBD). Prior digital analyses of normal breast tissues revealed that epithelial nuclear density (END) and TDLU involution are inversely correlated. Accordingly, we examined associations of END, TDLU involution, and MD in BBD clinical biopsies. This study included digitized images of 262 representative image-guided hematoxylin and eosin-stained biopsies from 224 women diagnosed with BBD, enrolled within the cross-sectional BREAST-Stamp project that were visually assessed for TDLU involution (TDLU count/100 mm2, median TDLU span and median acini count per TDLU). A digital algorithm estimated nuclei count per unit epithelial area, or END. Single X-ray absorptiometry of prebiopsy ipsilateral craniocaudal digital mammograms measured global and localized MD surrounding the biopsy region. Adjusted ordinal logistic regression models assessed relationships between tertiles of TDLU and END measures. Analysis of covariance examined mean differences in MD across END tertiles. TDLU measures were positively associated with increasing END tertiles [TDLU count/100 mm2, ORT3vsT1: 3.42, 95% confidence interval (CI), 1.87-6.28; acini count/TDLUT3vsT1, OR: 2.40, 95% CI, 1.39-4.15]. END was significantly associated with localized, but not, global MD. Relationships were most apparent among patients with nonproliferative BBD. These findings suggest that quantitative END reflects different but complementary information to the histologic information captured by visual TDLU and radiologic MD measures and merits continued evaluation in assessing cellularity of breast parenchyma to understand the etiology of BBD.

Relationship of circulating insulin-like growth factor-I and binding proteins 1-7 with mammographic density among women undergoing image-guided diagnostic breast biopsy.

BACKGROUND: Mammographic density (MD) is a strong breast cancer risk factor that reflects fibroglandular and adipose tissue composition, but its biologic underpinnings are poorly understood. Insulin-like growth factor binding proteins (IGFBPs) are markers that may be associated with MD given their hypothesized role in breast carcinogenesis. IGFBPs sequester IGF-I, limiting its bioavailability. Prior studies have found positive associations between circulating IGF-I and the IGF-I:IGFBP-3 ratio and breast cancer risk. We evaluated the associations of IGF-I, IGFBP-3, and six other IGFBPs with MD. METHODS: Serum IGF measures were quantified in 296 women, ages 40-65, undergoing diagnostic image-guided breast biopsy. Volumetric density measures (MD-V) were assessed in pre-biopsy digital mammograms using single X-ray absorptiometry. Area density measures (MD-A) were estimated by computer-assisted thresholding software. Age, body mass index (BMI), and BMI2-adjusted linear regression models were used to examine associations of serum IGF measures with MD. Effect modification by BMI was also assessed. RESULTS: IGF-I and IGFBP-3 were not strongly associated with MD after BMI adjustment. In multivariable analyses among premenopausal women, IGFBP-2 was positively associated with both percent MD-V (ß = 1.49, p value = 0.02) and MD-A (ß = 1.55, p value = 0.05). Among postmenopausal women, positive relationships between IGFBP-2 and percent MD-V (ß = 2.04, p = 0.003) were observed; the positive associations between IGFBP-2 and percent MD-V were stronger among lean women (BMI < 25 kg/m2) (ß = 5.32, p = 0.0002; p interaction = 0.0003). CONCLUSIONS: In this comprehensive study of IGFBPs and MD, we observed a novel positive association between IGFBP-2 and MD, particularly among women with lower BMI. In concert with in vitro studies suggesting a dual role of IGFBP-2 on breast tissue, promoting cell proliferation as well as inhibiting tumorigenesis, our findings suggest that further studies assessing the role of IGFBP-2 in breast tissue composition, in addition to IGF-1 and IGFBP-3, are warranted.

Using deep convolutional neural networks to identify and classify tumor-associated stroma in diagnostic breast biopsies.

The breast stromal microenvironment is a pivotal factor in breast cancer development, growth and metastases. Although pathologists often detect morphologic changes in stroma by light microscopy, visual classification of such changes is subjective and non-quantitative, limiting its diagnostic utility. To gain insights into stromal changes associated with breast cancer, we applied automated machine learning techniques to digital images of 2387 hematoxylin and eosin stained tissue sections of benign and malignant image-guided breast biopsies performed to investigate mammographic abnormalities among 882 patients, ages 40-65 years, that were enrolled in the Breast Radiology Evaluation and Study of Tissues (BREAST) Stamp Project. Using deep convolutional neural networks, we trained an algorithm to discriminate between stroma surrounding invasive cancer and stroma from benign biopsies. In test sets (928 whole-slide images from 330 patients), this algorithm could distinguish biopsies diagnosed as invasive cancer from benign biopsies solely based on the stromal characteristics (area under the receiver operator characteristics curve = 0.962). Furthermore, without being trained specifically using ductal carcinoma in situ as an outcome, the algorithm detected tumor-associated stroma in greater amounts and at larger distances from grade 3 versus grade 1 ductal carcinoma in situ. Collectively, these results suggest that algorithms based on deep convolutional neural networks that evaluate only stroma may prove useful to classify breast biopsies and aid in understanding and evaluating the biology of breast lesions.